#include <math.h>

int initphi(float *phi, int columns, int rows)
{
  int k, l;
  for(k=0;k<rows;k++)
    for(l=0;l<columns;l++)
      phi[k*columns+l] = 
	- sqrtf((columns/2-l)*(columns/2-l)
	      +(rows/2-k)*(rows/2-k))
	+ (columns<rows?columns:rows)/3;
  return 0;
}

